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Abstract - The Mars Microprobe mission will 
provide the first opportunity for subsurface 
measurements, including water detection, near 
the south pole of Mars. In this paper, per- 
formance of the Microprobe aeroshell design is 
evaluated through development of a six-degree- 
of- freedom (6-DOF) aerodynamic database and 
flight dynamics simulation. Numerous mission 
uncertainties are quantified and a Monte- Carlo 
analysis is performed to statistically assess mis- 
sion performance. Results from this 6-DOF 
Monte- Carlo simulation demonstrate that, in a 
majority of the cases (approximately 2-cr), the 
penetrator impact conditions are within current 
design tolerances. Several trajectories are iden- 
tified in which the current set of impact re- 
quirements are not satisfied. From these cases, 
critical design parameters are highlighted and 
additional system requirements are suggested. 
In particular, a relatively large angle-of-attack 
range near peak heating is identified. 
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1 Introduction 

The objective of NASA’s New Millennium pro- 
gram is to demonstrate and flight qualify tech- 
nology elements required for the science mis- 
sions of the next century [1]. The program’s 
second flight project, Deep Space Two (DS- 
2) is focused on the design of two small Mars 
entry probes. As a result, DS-2 is often re- 
ferred to as the Mars Microprobe mission. This 
proof-of-concept system is intended to demon- 
strate key elements of future network science 
missions [2, 3]. Attached to the cruise stage 
of the Mars 98 Surveyor Lander, these two 
Microprobe vehicles will be launched to Mars 
by a Delta II rocket in January 1999, arriving 
in December 1999. Each of these Microprobe 
capsules houses instrumented penetration de- 
vices designed to analyze the subsurface layers 
by performing soil sampling and water detec- 
tion. On impact, the penetrators are designed 
to pierce their protective aeroshells, driving this 
subsurface instrumentation 0. 3-2.0 m below the 
surface. Subsurface data will be relayed back to 
Earth through a link with the Mars Global Sur- 
veyor orbiter (September 1997 Mars arrival). 

The entry, descent, and impact (EDI) phase of 
the DS-2 mission begins as the two capsules are 
mechanically separated from the cruise stage [4]. 
This event is preceded by separation of the Mars 
98 Lander from the cruise stage (approximately 
1.5 s earlier). As a result of (1) the brief period 
between these two separation events, (2) the 
lack of control of the cruise-stage after the 98 
Lander separation, and (3) geometric mounting 



constraints which do not allow the Microprobe 
vehicles to be aligned with the flight path, the 
capsules will separate in an unknown angular 
orientation with non-zero angular rates. Sta- 
ble flight of the Microprobe vehicles must be 
achieved passively, and maintained until surface 
impact. 

Design of the DS-2 entry probes is compli- 
cated by several unique aerodynamic challenges. 
The vehicles must possess enough aerodynamic 
stability to achieve passive re- orient at ion from 
an arbitrary initial motion prior to peak heat- 
ing. Since stable flight at impact is required, 
the supersonic and transonic dynamic stability 
problems which have plagued other entry mis- 
sions [5, 6, 7] must also be mitigated. Addi- 
tionally, the penetrators must be protected from 
the intense aerothermodynamic environment of 
a 7.0 km/s Mars entry and satisfy a stringent 
set of surface impact constraints. 

In this paper, the criteria used to select the 
aeroshell geometry are presented. After re- 
view of the aeroshell shape and mass prop- 
erties, compilation of the Microprobe aerody- 
namic database is discussed. This database is 
compiled from past studies, computational fluid 
dynamic calculations, and ground-based test 
data. Development of a six-degree-of-freedom 
(6-DOF) Monte-Carlo trajectory simulation for 
Microprobe EDI is also presented. Results from 
this 6-DOF Monte-Carlo simulation are used to 
statistically assess the effect of combinatorial 
variations in the significant EDI parameters. 

2 Nomenclature 

A reference aerodynamic surface area, m 2 
b yaw /roll reference length, m 

c pitch reference length, m 

C aerodynamic force or moment coefficient 
g Mars surface slope, deg 

Kn Knudsen number 

m mass, kg 

M Mach number 

I body-axis moments of inertia, kg-m 2 

q stagnation-point heat rate, W/cm 2 


Q dynamic pressure, N/m 2 

V velocity, m/s 

a angle-of-attack, deg 

/3 ballistic coefficient kg/m 2 

7 flight-path or incidence angle, deg 

Subscripts 

a relative to the atmosphere 

A axial force 

D drag force 

1 static rolling moment 

m static pitching moment 

mq dynamic pitching moment 
n static yawing moment 

N normal force 

nr dynamic yawing moment 
p penetrator 

r relative to the horizon 

t total 

Y side force 

3 Surface Impact Requirements 

The penetrators have been designed to oper- 
ate properly under a range of impact condi- 
tions. Mission success demands that the EDI 
system meet several surface impact constraints. 
Three-cr requirements on surface impact veloc- 
ity (140 < V r < 200 m/s), penetration angle- 
of-attack (0 < a p < 10 degrees), and pen- 
etration incidence angle (| 7 P | < 20 degrees) 
have been specified [4]. These requirements 
are currently being validated through a rigor- 
ous ground-testing program. 

The dynamics of the surface impact event are 
illustrated in Fig. 1. The Microprobe aeroshell, 
local horizon, and local surface slope are shown, 
along with the velocities with respect to the 
ground ( V r ) and atmosphere (V a ). By conven- 
tion, the flight-path angles shown are negative. 
Atmospheric winds cause the difference between 
V r and V a . Total angle-of-attack (at in Fig. 1) 
is defined as the angle between the vehicle’s axis 
of symmetry and V a . 

During flight, the forces on the aeroshell are a 
function of the relationship between the vehicle 
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Figure 1 . Definition of Mars Microprobe sur- 
face impact angles. 

and the atmospheric velocity vector (at and V a )- 
However, at impact, the orientation of the Mi- 
croprobe payload relative to the surface is of sig- 
nificance. The penetration angle-of-attack and 
incidence angle are defined as: 

Otp = a t + ( 7 a - 7 r) (1) 

7p = 7 r + 90° - g (2) 

As discussed below, the surface impact velocity 
and incidence angle constraints may be achieved 
through the selection of the appropriate ballis- 
tic coefficient (see Section 5); whereas, satis- 
faction of the impact angle-of-attack constraint 
is a function of the vehicle’s aerodynamic sta- 
bility (geometry and center-of-gravity location). 
Aerodynamic design of the Microprobe capsules 
is presented in Section 4. The variance in each 
of the significant impact parameters as well as 
statistical data regarding the aeroshell heating 
environment is presented in Section 5. 

4 Analysis 

Aeroshell Selection 

Selection of the aeroshell for Mars Microprobe 
requires consideration of the unique objectives 
of the mission. A passive enclosure is required 
to safely deliver the penetrator payload through 
entry to impact with the surface. The aeroshell 
must decelerate the vehicle during its descent 
to a prescribed impact velocity. It must possess 
sufficient stability to correct any initial tum- 
bling motion to forward- facing flight early in the 
trajectory and maintain that orientation within 


small tolerances until impact. Finally, it must 
protect the payload from intense aerodynamic 
heating. To meet these objectives a 45-degree 
half-angle cone with rounded nose and shoul- 
ders is selected for the forebody. The afterbody 
is hemispherical with its center at the vehicle’s 
center-of-gravity location. 

Blunted 45-degree sphere-cones were used for 
the successful Pioneer- Venus and Galileo mis- 
sions [8]. Both of these missions entered at- 
mospheres much denser than Mars. For the 
Mars entries of Viking and Mars Pathfinder, 
70-degree sphere-cones with a zero angle-of- 
attack drag coefficient near 1.7 (versus the 45- 
degree cone value of 1.05) were selected [9]. 
Choice of cone angle calls for a compromise of 
drag, stability and packaging. Blunter cones ex- 
hibit more drag per surface area; sharper cones 
possess more stability. Viking and Pathfinder 
make use of high drag aeroshells since both 
of these entries required deceleration of much 
heavier spacecraft at sufficiently high altitudes 
for parachute deployment. In contrast, the 
Mars Microprobe vehicles are more than two 
orders of magnitude lighter and must impact 
the surface at a high velocity (140-200 m/s). 
Additionally, each Microprobe capsule requires 
the highest possible aerodynamic stability to re- 
cover quickly from any initial tumbling motion. 

The degree of nose bluntness has little effect 
on the drag coefficient for a 45-degree half- 
angle cone, although increased bluntness does 
slightly decrease static stability. On the other 
hand, increased bluntness decreases the stagna- 
tion point heat rate during the hypersonic por- 
tion of the trajectory. Selecting the appropri- 
ate degree of nose bluntness is a compromise 
of these factors. For Microprobe, a nose radius 
equal to half of the vehicle’s overall base radius 
is an acceptable value. This is the same ra- 
tio used in the Pioneer- Venus and Galileo entry 
probes [7, 10, 11]. Similarly, rounding the ve- 
hicle’s shoulders is performed to decrease local 
heating. Rounding the shoulders decreases both 
drag and stability. Again, the Pioneer- Venus 
value of shoulder radius equal to one tenth the 



nose radius is specified for Microprobe. Al- 
though it is possible to optimize the amount 
of nose and shoulder rounding for the specific 
Microprobe mission, the selection of previously 
used ratios appears adequate and also allows 
the use of an extensive body of existing aero- 
dynamic test and flight data. 

Selection of the hemispherical afterbody is 
based on the Planetary Atmosphere Experi- 
ments Test (PAET) probe [12]. The hemispher- 
ical afterbody specified for Microprobe serves 
two purposes. First, since the vehicle may 
be tumbling initially, it may encounter the at- 
mosphere while traveling backwards. A hemi- 
spherical afterbody with center at the vehicle’s 
center-of-gravity is not stable in this orienta- 
tion and will foster rotation to a forward fac- 
ing attitude. Second, this afterbody has been 
shown to decrease the dynamic instability ob- 
served in blunt vehicles traversing the transonic 
flight regime [13]. Regarding backwards stabil- 
ity, it is of interest to note that Pioneer- Venus, 
Galileo, Viking and Mars Pathfinder were all 
hypersonically stable in a backwards orienta- 
tion. To prevent this occurrence, each entry ve- 
hicle was oriented nose-first and spin-stabilized 
to assure a forward- facing attitude at the at- 
mospheric interface. Spin stabilization is not 
an option for Microprobe; however, the hemi- 
spherical afterbody assures the vehicle will not 
trim in a backwards-facing attitude. 

The geometry of the Mars Microprobe aeroshell 
is depicted in Fig. 2. As shown, a 45-degree 
sphere cone with nose radius of 0.0875 m, shoul- 
der radius of 0.00875 m, and maximum radius 
of 0.175 m has been selected. The afterbody 
shape is a hemispherical section with radius of 
0.183 m centered about the vehicle’s center-of- 
gravity. The center-of-gravity is located 0.0902 
m aft of the nose on the vehicle’s symmetry axis. 

Aerodynamics 

The Mars Microprobe aerodynamic database 
was derived from a combination of computa- 
tional fluid dynamic (CFD) calculations and 



Figure 2. Mars Microprobe aeroshell geometry 
and mass properties. 
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Figure 3. Sources used to assemble Mars Mi- 
croprobe aerodynamic database. 

ground-based test data. A detailed description 
of the vehicle’s aerodynamic characteristics is 
provided in Ref. [14]. Sources of the static aero- 
dynamic predictions are illustrated in Fig. 3. 
Free molecular and Direct Simulation Monte 
Carlo (DSMC) computations were performed to 
characterize the rarefied and transitional flow 
regimes. These results were supplemented by 
thermochemical nonequilibrium computational 
fluid dynamic calculations obtained with the 
Langley Aerothermodynamic Upwind Relax- 
ation Algorithm (LAURA in Fig. 3) in the con- 
tinuum hypersonic flow regime. This analysis 
tool was extensively used in the prediction of the 
Mars Pathfinder aerodynamics [15]. Pioneer- 
Venus wind tunnel data was used in the super- 
sonic, transonic, and subsonic regimes [16, 17]. 



Validation and extrapolation of these exist- 
ing results was made possible through addi- 
tional computational solutions obtained in the 
transonic and subsonic flight regimes with the 
Thin-Layer Navier-Stokes 3-Dimensional pro- 
gram [18] (TLNS3D in Fig. 3). 

Dynamic damping coefficients were extracted 
from Pioneer- Venus and Viking wind tunnel 
test data [5, 6, 7]. In addition, transonic bal- 
listic range data was produced on a Micro- 
probe model with the correct center-of-gravity 
location [19, 20]. Dynamic stability estimates 
could not be obtained computationally within 
the time constraints of the present analysis. Be- 
cause of the transonic dynamic instability prob- 
lems which have plagued other entry vehicle de- 
signs, additional data is being gathered in a 
pressurized facility in which the flight Reynolds 
number can be duplicated [19]. 

To produce a cohesive database from these di- 
verse sources, modification of the original data 
set was required [21]. A bridging function 
(shaped by the DSMC results) was used in the 
transitional region between the free molecular 
and continuum results. Explicit calculation of 
the transitional aerodynamics by DSMC meth- 
ods, although possible, is computationally pro- 
hibitive. Instead, selected DSMC results were 
used to anchor and shape the bridging function. 
This function provides a smooth variation of the 
vehicle’s aerodynamic characteristics based on 
the free molecular and continuum hypersonic 
computations. 

Within the database [21], the continuum hy- 
personic aerodynamics were assumed to vary 
with angle of attack in a similar manner to 
that predicted by Newtonian flow. This New- 
tonian variation was then scaled to reproduce 
specific LAURA computations obtained at 0 
and 10-degrees angle of attack. The supersonic 
and transonic wind tunnel data overlapped, so 
the two sets were blended in the Mach number 
region between 1.65 and 2.16. Similarly, the 
Pioneer- Venus and Viking dynamics data were 
blended in the Mach number region between 1.2 
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Figure 4. Mars Microprobe 0-degree angle-of- 
attack drag coefficient. 

and 1.3. Finally, to account for small differences 
in the aeroshell geometry, the axial force wind 
tunnel values were scaled to match the compu- 
tational results. 

Early versions of the aerodynamics routine used 
simple linear interpolation, providing value- 
continuity between segments. Although a 
twice-differentiable database was sought, pro- 
viding this level of continuity at the data 
points resulted in unacceptable behavior be- 
tween the data. As a compromise, an overlap- 
ping parabola technique, which provides slope- 
continuity, was used. As the FORTRAN rou- 
tine was developed, care was taken to minimize 
memory overhead. Furthermore, since the rou- 
tine is called many times by 6-DOF POST, an 
effort was made to create a computationally ef- 
ficient algorithm. For a given flight condition 
and vehicle angular orientation, the database 
provides estimates of C U, Cy, Cy, C m , C n , C/, 
C m<7 , and C nr for use in the 6-DOF trajectory 
simulation [21]. 

The 0-degree angle-of-attack Microprobe drag 
coefficient is shown in Fig. 4 as a function 
of Knudsen (Kn) and Mach (M) numbers. 
Knudsen number is defined as the ratio of the 
gas’ mean free path to the vehicle’s diameter. 
This similarity parameter is used as the in- 
dependent variable in the rarefied and transi- 




tional aerodynamic regimes. Initially, Kn will 
be large. For values larger than 10, the aero- 
dynamic forces are computed solely from the 
free molecular flow solutions. Free molecular 
flow assumes there are no collisions between gas 
molecules in the flow field. Unlike hypersonic 
continuum aerodynamics (where forces exerted 
on the body are essentially the integrated ef- 
fect of surface pressures alone), free molecular 
flow aerodynamics contain a significant shear 
stress contribution. As the entry proceeds into 
the upper atmosphere, both the mean free path 
and Knudsen number decrease and collisions 
between particles must be taken into account. 
In this regime, where 0.001 < Kn < 10, the 
aerodynamics are computed from the DSMC- 
anchored bridging function. 

As lower altitudes are reached (below 55 km al- 
titude for Microprobe), the Knudsen number 
drops below 0.001 and the continuum meth- 
ods are used to compute vehicle aerodynamics. 
Here, Mach number is the appropriate aerody- 
namic similarity parameter. Fig. 4 shows that 
drag coefficient (at a given angle of attack) is 
approximately constant above Mach 5. The in- 
crease in drag coefficient supersonically, is a re- 
sult of the sonic line shifting from the nose re- 
gion to the shoulder region of the aeroshell be- 
tween Mach 5 and 2. This shift has a signifi- 
cant impact on the pressure distribution caus- 
ing axial force (equivalent to drag at 0-degree 
angle-of-attack) to increase while normal force 
and moment coefficient decrease. Transonically, 
drag coefficient decreases. Here the data of 
Refs. [17, 20] and the TLNS3D computational 
solutions are used. The axial force data of 
Ref. [17] is decreased to account for the Micro- 
probe configuration’s large hemispherical after- 
body. 

Atmospheric Flight Dynamics 

Six- degree- of- freedom (6-DOF) trajectory anal- 
ysis is performed using the Program to Opti- 
mize Simulated Trajectories (POST) [22]. This 
program has been used previously in the devel- 
opment of the Mars Pathfinder entry, descent, 


and landing (EDL) strategy [23, 24]. Six-DOF 
POST is also being used by the Pathfinder EDL 
operations team [25, 26]. In the present study, 
POST is used to numerically integrate the 6- 
DOF equations of motion from a given entry 
state to surface impact. An eighth-order Runge- 
Kutta integration technique is employed [27]. 
The Microprobe aerodynamic database as well 
as Mars atmospheric, gravitational, and surface 
models are inputs to the simulation. Atmo- 
spheric modeling for this mission is hampered 
by the lack of surface measurements for the 
southern hemisphere of Mars (the target im- 
pact site is 73-77 degrees South latitude). Using 
projected data obtained from the Viking lan- 
ders and a global circulation model, Zurek and 
Richardson have constructed nominal and per- 
turbed atmospheric profiles [28]. These models 
are used in the current simulation. 

In the present analysis, uncertainties are applied 
in all simulation model inputs. These uncertain- 
ties arise from numerous sources including (1) 
technology limitations (e.g., current interplane- 
tary navigation or mass-balance accuracies), (2) 
a lack of knowledge concerning the Mars atmo- 
sphere, (3) computational or measurement un- 
certainty associated with the aerodynamic anal- 
yses, and (4) unknown separation orientation 
and angular rate. Therefore, in this analysis, 
an attempt was made to quantify and model 
the degree of uncertainty in each of 29 major 
parameters. 

The uncertainty range attributed to each of 
these parameters is listed in Table 1. For a 
parameter with more than one variance (e.g., 
aerodynamics or winds), the uncertainty is es- 
timated using linear interpolation between the 
regions given in Table 1. Gaussian distributions 
are sampled for most parameters. However, the 
initial orientation, center-of-gravity offset quad- 
rant, and wind direction quadrant are deter- 
mined from uniform distributions. The topog- 
raphy variation is modeled by a non-symmetric 
Gaussian distribution centered at 5 km. This 
distribution is illustrated in Fig. 5. 



Table 1. 6-DOF Monte-Carlo variables. 


i oo r 


EDI Parameter 

Nominal 

Value 

3-cr 

Variance 

Initial state, 7 , deg 

-13.25 

± 0.4 

Initial at, deg 

90.0 

=b 90.0 

Initial pitch rate, deg/s 

6.0 

zb 5.0 

Initial yaw rate, deg/s 

6.0 

=b 5.0 

Initial roll rate, deg/s 

6.0 

=b 5.0 

X-axis eg position, mm 

90.2 

zb 5.0 

X-axis eg offset, mm 

0.0 

=b 1.0 

Mass, kg 

2.73 

=b 0.273 

l xx , kg-m 2 

0.0106 

zb 0.0003 

lyy, Izz, kg-m 2 

0.0105 

=b 0.0003 

Ut /5 Icczj Il/Zj kg-m 

0.0 

=b 0.0003 

Ca, Kn > 0.1 

See Fig. 4 

± 10 % 

C N, Kn > 0.1 

See Ref. [14] 

=b 0.10 

Cm, Kn > 0.1 

See Ref. [14] 

=b 0.006 

Ca, Kn < 0.1, M > 10 

See Fig. 4 

± 2% 

Cat, Kn < 0.1, M > 10 

See Ref. [14] 

zb 0.05 

Cm, Kn < 0 . 1 , M > 10 

See Ref. [14] 

=b 0.003 

Ca, Kn < 0.1, M < 5 

See Fig. 4 

± 10 % 

Cat, Kn < 0.1, M < 5 

See Ref. [14] 

=b 0.10 

Cm, Kn < 0 . 1 , M < 5 

See Ref. [14] 

=b 0.006 

C mq and Cnr, M > 6 

See Ref. [14] 

zb 20% 

C mq and Cnr, M < 3 

See Ref. [14] 

-50, +110% 

Density 

See Ref. [28] 

See Fig. 6 

Temperature 

See Ref. [28] 

± 150% 

Wind above 50 km 

See Ref. [28] 

See Fig. 7 

Wind below 10 km 

See Ref. [28] 

See Fig. 7 

Wind gust, m/s 

0.0 

± 30.0 

Surface altitude, km 

5.0 

-4, +1 

Surface slope, deg 

0.0 

± 5.0 



Figure 5. Impact altitude distribution. 



Figure 6. 3 -a density variation used in the 

6-DOF Monte-Carlo simulation. 



Wind speed, m/s 


Figure 7. 3 -a wind profiles used in the 6-DOF 
Monte-Carlo simulation. 

The 3 -a variances in atmospheric density and 
winds are shown in Figs. 6 and 7. As illustrated 
in these figures, three atmospheric profiles with 
varying levels of visible column depth, or opac- 
ity (0.05, 0.2, 0.5), are used in the present sim- 
ulation. Within the Monte-Carlo simulation, 
these three nominal and perturbed atmospheric 
profiles are sampled evenly. 

5 Results And Discussion 

Impact Sizing and Nominal Trajectory 

For zero angle-of-attack flight of the Microprobe 
capsules, the impact conditions are completely 
determined by ballistic coefficient (/?) and sur- 





Figure 8. Mars Microprobe ballistic coefficient 
requirements, surface altitude = 5 km. 


18 

17 

16 

15 

> 14 

deg 

13 

12 

11 

10 



probe impact conditions, f3 = 27 kg/m 2 . 


face altitude. Assuming a surface slope of 5 de- 
grees, an altitude of 5 km, the mean density 
profile [28], and negligible wind speed, the fol- 
lowing values are determined from Fig. 8. For 
an impact speed of 160 m/s, a ballistic coeffi- 
cient of 27 kg/m 2 is required. Furthermore, the 
ballistic coefficient must be within the range of 
20.6 < f3 < 47 kg/m 2 to satisfy the impact ve- 
locity criterion (140 < V r < 200 m/s), and be- 
low 33.1 kg/m 2 to satisfy the penetration in- 
cidence angle constraint (| 7 P | < 20 degrees). 
Hence, the Microprobe design space is restricted 
such that 20.6 < (3 < 33.1 kg/m 2 . 

The design approach selected by the Mars Mi- 
croprobe project office is to baseline the largest 
diameter aeroshell that does not adversely im- 
pact the Mars 98 Lander. This approach, which 
allows for modest mass growth, yields the 350 
mm maximum diameter depicted in Fig. 2. As 
discussed earlier, a 45-degree sphere-cone has a 
continuum hypersonic Co of 1.05. Thus, for the 
nominal value of f3 = 27 kg/m 2 , a Microprobe 
mass of 2.73 kg is required. Without ballast, the 
current system mass estimate is 2.63 kg. This 
yields a mass margin of 3.8% at the nominal 
impact speed. As mentioned above, the impact 
constraints can be met with a ballistic coeffi- 
cient as high as (3 = 33.1 kg/m 2 ; hence, the Mi- 
croprobe heatshield is being sized for this value. 
At this upper limit of (3, the current design mass 
margin is 27.1%. 


The impact of surface altitude on the Fig. 8 
trade-space is shown in Fig. 9. For the impact 
site of interest, a surface altitude of 1 to 6 km 
above the Mars reference ellipsoid may be en- 
countered [4]; however, a majority of this terrain 
is thought to be 5 km above the Mars reference 
ellipsoid (see Fig. 5). As shown in Fig. 9, im- 
pact velocity changes approximately 5 m/s for 
each km of surface altitude and the penetration 
incidence angle decreases with surface altitude 
(at a rate of about 1.4 deg/km). For (3 = 27 
kg/m 2 , the impact velocity criterion is satisfied 
above 1.6 km altitudes; whereas the penetra- 
tion incidence angle constraint is satisfied over 
the complete range of expected surface altitude. 

The nominal Microprobe EDI trajectory is 
shown in Fig. 10. As shown, the atmospheric 
entry velocity is 6.9 km/s. The atmospheric 
interface is defined at a radius of 3522.2 km 
(surface altitude of 141.8 km). The peak de- 
celeration of 12.6 Earth g’s is achieved at an 
altitude of 44.1 km. Surface impact occurs 
roughly 290 sec after atmospheric interface, at 
160 m/s (Mach 0.66). Fig. 11 shows that the 
peak stagnation- point heat rate of 160 W/cm 2 
is achieved at approximately 80 sec (altitude of 
49.2 km) and is followed by peak dynamic pres- 
sure (3380 N/m 2 ) at approximately 94 sec into 
the EDI sequence. At impact, the integrated 
stagnation- point heat load is 7165 J/cm 2 . 
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Figure 10. Mars Microprobe nominal trajec- 
tory. 
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Figure 11 . Stagnation-point heat-rate and dy- 
namic pressure profiles. 

The nominal entry shown in Figs. 10 and 11 as- 
sumes a zero angle-of-attack entry interface ori- 
entation, mean atmosphere without wind, and 
no error in the center-of-gravity position, mass 
properties, or aerodynamic modeling. Figure 12 
presents the total angle-of-attack history for a 
similar entry initiated with a 90-degree angle- 
of-attack. As shown, total angle-of-attack is 
damped below 10 degrees within the first 50 
sec of atmospheric flight (by 76 km altitude) 
and continues to decrease as the dynamic pres- 
sure (Q) builds (see Fig. 11). Thus, although 
high angles-of-attack may occur early in the 
mission, the Microprobe aeroshells possess suf- 
ficient aerodynamic stability to provide passive 
re-orientation while in the upper atmosphere. 



0 50 100 150 200 250 300 


Time, s 

Figure 12. Mars Microprobe total angle-of- 
attack profile, initial 90-degree angle-of-attack. 

This hypersonic re- orient at ion issue is discussed 
further in the final section of this paper. 

In addition to demonstrating the hypersonic 
re-orientation capability of the Microprobe 
aeroshell, Fig. 12 also demonstrates that stat- 
ically stable flight can be maintained passively 
throughout EDI. In fact, for this case, the angle- 
of-attack at impact (at) is 1.5 degrees. Note 
that the presence of the transonic dynamic in- 
stability is evident in Figure 12. This phenom- 
ena is the cause for the small increase in total 
angle-of-attack from 225 to 275 sec. 

Because of Microprobe’s small moments of in- 
ertia (see Fig. 2), large angular rates are likely 
during the entry. This high frequency motion is 
evident in Fig. 12, particularly when dynamic 
pressure is large. Pitch and yaw rates as high 
as 200 deg/sec and an angle-of-attack oscilla- 
tion frequency of approximately 5 Hz are likely 
near peak Q. In comparison, Mars Pathfinder’s 
angular motion is characterized by pitch and 
yaw rates two orders of magnitude lower in 
the peak dynamic pressure region and a much 
lower angle-of-attack frequency [23]. Although 
Microprobe’s angular frequency diminishes as 
dynamic pressure decreases, an increased fre- 
quency occurs during transonic flight as a result 
of the vehicle’s dynamic instability. At impact, 
the vehicle’s pitch and yaw rates are more than 
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Figure 13. Impact velocity distribution. 


an order of magnitude lower than at peak dy- 
namic pressure and the angle-of-attack oscilla- 
tion frequency is approximately 2.5 Hz. 

Monte- Carlo Simulation 

During flight of the Mars Microprobe space- 
craft, a combination of off-nominal effects is 
likely to be encountered. Hence, it is impor- 
tant to statistically assess the effect of combi- 
natorial variations in all of the EDI parame- 
ters listed in Table 1. To accomplish this, two- 
thousand off-nominal cases were randomly esti- 
mated and simulated in a Monte- Carlo fashion. 
A 99.7% probability exists that each random pa- 
rameter will remain within the 3-cr uncertainty 
bounds of Table 1. In addition to a detailed set 
of impact conditions, the total angle-of-attack 
was monitored at discrete points along the heat 
pulse. Peak deceleration, stagnation- point heat- 
rate, and integrated heat load were also moni- 
tored. 

Impact Conditions-JUstogidims of the probable 
ranges in impact velocity, penetration angle-of- 
attack and surface incidence angles are shown 
in Figs. 13 - 15. Monte-Carlo statistics for 
these impact parameters are tabulated in Ta- 
ble 2. These data indicate that the current set 
of Microprobe impact requirements are not sat- 
isfied in a 3-cr sense. In particular, the 3-cr low 
impact velocity criteria (140 m/s) and the 3- 
cr high penetration incidence angle (20 degrees) 
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Figure 14. Penetration angle-of-attack distri- 
bution. 
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Figure 15. Penetration incidence angle distri- 
bution. 


Table 2. 6-DOF Mars Microprobe Monte- 
Carlo analysis results. 


EDI Parameter 

Mean 

Value 

3-cr 

Variance 

Impact velocity, V r , m/s 

156.1 

± 22.0 

Pen. angle-of-attack, a p , deg 

4.6 

+ 14.7 

Pen. incidence angle, j p , deg 

8.7 

zb 14.4 

Impact Mach number, km 

0.65 

=b 0.09 

Impact time, sec 

298.0 

=b 29.8 

Impact latitude, deg 

-74.7 

zb 1.7 

Impact longitude, deg 

148.0 

zb 1.1 

Peak q, W/cm 2 

154.1 

zb 11.9 

Peak deceleration, Earth g’s 

12.4 

=b 1.6 

Integrated heat load, J/cm 2 

6875. 

=b 416.7 

Peak q angle-of-attack, at, deg 

13.1 

+ 27.0 

Peak Q angle-of-attack, at, deg 

7.8 

+ 4.2 





Figure 16. Penetration angles dispersion. 
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Figure 17. Impact Mach number and velocity 
dispersion. 




Figure 18. Impact altitude and time disper- 
sion. 
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Figure 19. Impact footprint. 


are mildly violated. Note that these two crite- 
ria are currently satisfied to a 2 -a probability 
level. While a small increase in ballistic coeffi- 
cient could be used to adjust the impact velocity 
range, such an increase is not recommended as 
this would only exacerbate the penetration in- 
cidence angle problem (see Fig. 8). 

Distributions of the significant impact parame- 
ters are shown in Figs. 16 - 19. The two pen- 
etration angles show a significant degree of clus- 
tering for a p < 4.0 degrees and 2.0 < < 13.0. 

Note that the variation exhibited in Fig. 17 be- 
tween impact velocity and Mach number is a 
result of atmospheric temperature variability. 

The Microprobe impact footprint is presented in 
Fig. 19. This footprint should bear some resem- 


blance to the Mars 98 Lander footprint, since 
both spacecraft approach Mars on the same in- 
terplanetary trajectory. As a result of Micro- 
probe’s lower ballistic coefficient, the range of 
impact sites shown in Fig. 19 does not extend as 
far southeast as the Mars 98 Lander footprint. 
However, Microprobe’s 180 x 20 km ellipse is 
similar in downrange and larger in crossrange 
(Mars 98 Lander is actively controlled) than the 
predicted Mars 98 Lander ellipse [29]. 

Peak stagnation- point heating, integrated heat 
load, and peak atmospheric deceleration statis- 
tics are also presented in Table 2. Histograms of 
these EDI parameters are presented in Figs. 20- 
22. The variation in intensity of this aerother- 
modynamic environment does not currently 
constrain the mission design space because the 




Figure 23. Penetration angles of attack and 
Figure 20. Peak stagnation- point heat rate incidence probability distribution, 
distribution. 
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Figure 21. Integrated heat load distribution. 



10.5 11.0 11.5 12.0 12.5 13.0 13.5 14.0 14.5 15.0 
Peak deceleration, Earth g's 

Figure 22. Peak atmospheric deceleration dis- 
tribution. 


heatshield is designed to withstand a peak heat 
rate of approximately 200 W/cm 2 and an in- 
tegrated heat load of more than 8550 J/cm 2 , 
based on entry with a (3 of 33.1 kg/m 2 . Sim- 
ilarly, the deceleration levels shown in Fig. 22 
during atmospheric flight pale in comparison to 
30000 Earth g’s expected at impact. 

As presented in Table 2, the mean angle-of- 
attack at peak heating and impact do not ap- 
pear too severe. Since total angle-of-attack is 
a one-sided function, these variables do not ex- 
hibit a Gaussian distribution. Hence, symmet- 
ric 3-cr variances are not listed for angle-of- 
attack parameters. Instead, probability curves 
have been generated and are presented in 
Fig. 23 for the two significant impact angles (a p 
and t p ). These curves give the probability that 
the angle-of-attack will be below a given value. 
Hence, from Figs. 23, a 93% probability (less 
than 2 -a) currently exists that the penetration 
angle-of-attack will be below 10 degrees. Simi- 
larly, there is a 96% probability (slightly more 
than 2 -a) that the penetration incidence angle 
will be below 20 degrees. Note that as a result 
of surface slope variability, the penetration inci- 
dence angle may be slightly negative (as much 
as -2.78 degrees was obtained). 

The variances in all of the impact conditions 
may be reduced through improved modeling of 
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Figure 24. Angle-of-attack probability dis- 
tribution at several points along the increasing 
side of heat pulse. 

the surface winds, altitude, and slope. For ex- 
ample, most of the impact velocity variance 
observed in Fig. 13 is a result of surface al- 
titude uncertainty. Similarly, a large percent- 
age of the penetration angle variances presented 
in Figs. 14 and 15 are the result of the ± 5- 
degree surface slope currently assumed. The ex- 
pected range in these penetration angles could 
also be reduced by moving the center-of-gravity 
forward. 

Hypersonic Re- orientation - Through this analy- 
sis, a relatively large angle-of-attack range near 
peak heating has also been identified. As yet, 
heating calculations at this flight condition have 
only been performed at angles of attack below 
10 degrees. However, in this analysis, angles- 
of-attack as high as 45 degrees were obtained 
in the peak heating region. In fact, the mean 
angle-of-attack at peak heating (see Table 2) is 
13.07 degrees. These results have necessitated 
the performance of additional afterbody heating 
analyses. 

To bound this aerothermodynamic assessment, 
angle-of-attack probability functions at discrete 
points along the increasing side of heat pulse 
were generated from the present set of Monte- 
Carlo results. This data is shown in Fig. 24. 
Since EDI is initiated from a random angular 
orientation (0 < at < 180 degrees), this fig- 



Figure 25. Statistical variation in angle of at- 
tack as a function of stagnation- point heat rate. 


ure also provides a measure of the aeroshell’s 
hypersonic re- orient at ion capability. As a re- 
sult of the aeroshell’s aerodynamic stability, this 
re- orient at ion occurs relatively early in the at- 
mospheric flight. For example, by the time a 
q of 10 W/cm 2 is achieved, the angle-of-attack 
range has converged such that there is a 50% 
chance that the angle-of-attack will be below 30 
degrees. On the nominal trajectory, this heat 
rate occurs just 28.2 sec past the atmospheric 
interface at an altitude of 99.1 km. 

By the time 40 W/cm 2 (45.1 sec, 78.8 km alti- 
tude) is achieved, there is a 70% chance that the 
angle-of-attack will be below 20 degrees and a 
96% chance that the angle-of-attack will be be- 
low 30 degrees. As 100 W/cm 2 is approached 
(60.6 sec, 62.8 km altitude), the probability is 
95% that the angle-of-attack will be below 20 
degrees. This information is shown as a func- 
tion of stagnation- point heat rate in Fig. 25. 
In this figure, the peak heating region is at 150 
W/cm 2 . At peak heating, there is a 99% proba- 
bility that the angle-of-attack will be below 30.0 
degrees and a 90% probability that the angle- 
of-attack will be below 18.0 degrees. 

The vehicle’s hypersonic re- orient at ion capabil- 
ity is largely a function the vehicle’s aerody- 
namic shape and center-of-gravity position. For 
a given vehicle configuration, the damping pro- 
vided by the atmospheric forces depends on dy- 




namic pressure, Q. Unfortunately, as shown 
in Fig. 11, peak Q occurs approximately 15 
sec past peak q. Hence, the forces providing 
this passive angle-of-attack convergence are not 
at their greatest level until beyond peak heat- 
ing. In fact, in 97% of the cases simulated, 
the angle-of-attack stays below 10-degrees from 
peak dynamic pressure to Mach 2 (see Fig. 12). 
Improvements in the hypersonic re- orient at ion 
level prior to peak q can be achieved by moving 
the center-of-gravity position forward, reducing 
the products of inertia, or reducing the angular 
rates associated with the initial separation ma- 
neuver. This is currently the number one issue 
confronting the EDI team. In an effort to de- 
termine the mission impact of this phenomenon, 
additional afterbody heating analyses are being 
performed. 

6 Conclusions 

Scheduled for flight in 1999, the New Millen- 
nium Mars Microprobe mission will provide the 
first opportunity for subsurface measurements, 
including water detection, near the south pole 
of Mars. Design of the Microprobe aeroshells 
poses several unique aerodynamic challenges in- 
cluding passive hypersonic re- orient at ion of an 
initially tumbling body and stringent stability 
constraints during a subsonic impact. The cri- 
teria used in the selection a 45-degree sphere- 
cone forebody with hemispherical afterbody are 
presented in this paper. 

The performance of the Microprobe aeroshell 
design is then evaluated through the develop- 
ment of a six-degree-of- freedom (6-DOF) aero- 
dynamic database and flight dynamics simula- 
tion. Numerous mission uncertainties are quan- 
tified and a Monte- Carlo analysis is performed 
to statistically assess mission performance. Re- 
sults from this 6-DOF Monte-Carlo simulation 
demonstrate that, in a majority of the cases 
(generally 2-cr), the penetrator impact condi- 
tions are within current design tolerances. In 
particular, the 3 -a range of impact velocity is 
134-178 m/s. The 2-cr high penetration angle- 
of-attack is 12 degrees and the 2-cr high pene- 


tration incidence angle is 18 degrees. The Mi- 
croprobe impact footprint extends 180x20 km 
with its center at -74.73 degrees South latitude, 
147.98 degrees East longitude. Suggestions for 
improvement are made to enhance Microprobe 
EDI performance. 

In addition, a relatively large angle-of-attack 
range near peak heating is identified which 
has necessitated the performance of additional 
afterbody heating analyses. These afterbody 
heating analyses are strongly coupled to the 
hypersonic re- orient at ion performance derived 
from the present results. In particular, a 50% 
probability exists that the angle-of-attack will 
be below 30 degrees before a q of 10 W/cm 2 
is achieved; whereas, a 96% probability exists 
that the angle-of-attack will be below 30 de- 
grees before q = 40 W/cm 2 . As 100 W/cm 2 
is approached, there is a 95% chance that the 
angle-of-attack will be below 20 degrees. At 
peak heating, a mean angle-of-attack of 13.07 
degrees was obtained and a 90% probability ex- 
ists that this parameter will be below 18 de- 
grees. 
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